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ABSTRACT 

rn" 

The distribution in color and absolute magnitude is presented for 55,158 galaxies 
CN ■ taken from the Sloan Digital Sky Survey in the redshift range 0.08 < z < 0.12, as a 

function of galaxy number overdensity in a line-of-sight cylinder of transverse radius 
1 h^ 1 Mpc. In all environments, bulge-dominated galaxies (defined to be those with 
C | radial profiles best fit with large Sersic indices) form a narrow, well defined color- 

magnitude relation. Although the most luminous galaxies reside preferentially in the 
highest density regions, there is only a barely detectable variation in the color (zero- 
■ point) or slope of the color-magnitude relation (< 0.02 mag in 0,1 [<? — r] or [B — V]). 

These results constrain variations with environmental density in the ages or metallicities 
of typical bulge-dominated galaxies to be under 20 percent. 

o . 

Subject headings: cosmology: observations — galaxies: clusters: general — galaxies: 
elliptical and lenticular, cD — galaxies: evolution — galaxies: fundamental parameters 
- galaxies: statistics 
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1. Introduction 

Red, bulge-dominated galaxies show many precise regularities, including narrow, mass-dependent 
distributions of colors, surface-brightnesses, radial profile shapes, ages, chemical abundances, ve- 
locity dispersions, and mass-to-light ratios (e.g., Baum 1959; Faber 1973; Faber & Jackson 1976; 
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Visvanathan & Sandage 1977; Djorgovski h Davis 1987; Dressier et al. 1987; Kormendy & Djor- 
govski 1989; Roberts k, Haynes 1994; Burstein et al. 1997; Terlevich et al. 2001; Eisenstein et al. 
2003; Bernardi et al. 2003a, b,c; Blanton et al. 2003a). The most massive of these bulge-dominated 
galaxies reside preferentially in the most dense environments, and especially in rich clusters of 
galaxies (Dressier 1980; Postman Sz Geller 1984); indeed the most massive red galaxy populations 
have the highest average environmental overdensities (Hogg et al. 2003). 

Are these trends, i.e., the relationships between mass and environment and between mass 
and everything else, related? In this Letter we look at the color-magnitude relationship of red, 
bulge-dominated galaxies as a function of environment. 

Generically, in models in which structure grows gravitationally, the most dense regions of the 
Universe will have collapsed earlier and will contain the most massive objects. Detailed modeling of 
galaxy formation in current cosmological simulations bears this out; galaxies in very high density 
regions are predicted to be more luminous and more red (older) than those in typical density 
regions (e.g., Cen & Ostriker 1993; Kauffmann et al. 1993; Kauffmann 1996; Kauffmann et al. 
1999; Benson et al. 2000; Diaferio et al. 2001). While these models predict morpholgy, color, and 
luminosity variations with overdensity, it is not clear whether they predict color differences at fixed 
morphology and luminosity. 

The Sloan Digital Sky Survey (SDSS) is the best available data set for work on these questions, 
because of its sample size, high signal-to-noise imaging, sky coverage, and complete spectroscopy 
(e.g., York et al. 2000). Indeed, the SDSS has already made some progress on related issues, 
including the dependence of clustering on luminosity and color (Zehavi et al. 2002), the star- 
formation rate as a function of environment (Gomez et al. 2003), the mean red-galaxy spectrum as 
a function of environment (Eisenstein et al. 2003) , and the mean galaxy environment as a function 
of color and luminosity (Hogg et al. 2003; Blanton et al. 2003a). 

In what follows, a cosmological world model with (Om,^a) = (0.3,0.7) is adopted, and the 
Hubble constant is parameterized Hq = 100 h kms _1 Mpc" 1 , for the purposes of calculating dis- 
tances and volumes (e.g., Hogg 1999). 

2. Data 

The SDSS is taking ugriz CCD imaging of 10 4 deg 2 of the Northern Galactic sky, and, from 
that imaging, selecting 10 6 targets for spectroscopy, most of them galaxies with r < 17.77 mag (e.g., 
Gunn et al. 1998; York et al. 2000; Stoughton et al. 2002). All the data processing: astrometry (Pier 
et al. 2003), source identification, deblending and photometry (Lupton et al. 2001), calibration 
(Fukugita et al. 1996; Smith et al. 2002), spectroscopic target selection (Eisenstein et al. 2001; 
Strauss et al. 2002; Richards et al. 2002), spectroscopic fiber placement (Blanton et al. 2003), and 
spectroscopic data reduction are performed with automated SDSS software. Redshifts are measured 
on the reduced spectra by an automated system, which models each galaxy spectrum as a linear 
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combination of stellar populations (Schlegel, in preparation). 

To each galaxy's extracted radial profile, a seeing-corrected Sersic model is fit (as in Blanton 
et al. 2003a). The Sersic model is that surface brightness I is proportional to exp(— r/ro) 1 /™, 
where r is the angular radius from the galaxy center, ro is a characteristic radius, and n is a free 
parameter known as the "Sersic index" . Galaxies with more concentrated radial profiles have larger 
n; an exponential disk has n = 1 and a de Vaucouleurs profile has n = 4. 

Galaxy luminosities (absolute magnitudes) and colors (measured by the standard SDSS pet- 
rosian technique; Petrosian 1976) are computed in fixed bandpasses, using Galactic extinction cor- 
rections (Schlegel et al. 1998) and K corrections (computed with kcorrect vl_ll; Blanton et al. 
2003). Galaxy magnitudes are K corrected not to the redshift z = observed bandpasses but to 
bluer bandpasses 0A g, 0A r and 01 i "made" by shifting the SDSS g, r, and i bandpasses to shorter 
wavelengths by a factor of 1.1 (as in Blanton et al. 2003, a). This means that galaxies at redshift 
z = 0.1 have trivial (but not zero; Hogg et al. 2002) K corrections. The full error analysis of SDSS 
photometry is not complete, but scatter in the g, r, and i band magnitudes appears to be better 
than 0.03 mag (RMS). 

For a red galaxy, 0,1 [g — r] (calibrated to the AB magnitude system) can be converted approx- 
imately to [B — V] with 

[B - V] AB = 0.83 ai [g - r] - [0.025 mag] 

[B~V] AB = [B-V] yega - [0.115 mag] , (1) 

where the subscripts indicate AB or Vega-relative calibration. This transformation is robust because 
the 01 g and 01 r bandpasses are very close to B and V in effective wavelength. 

The sample of galaxies used here is a subset of an SDSS statistical sample known as "NYU 
LSS samplel2", further selected to have apparent magnitude in the range 14.5 < r < 17.77 mag, 
and redshift in the narrow range 0.08 < z < 0.12. These cuts left 55,158 galaxies. 

The sample is statistically complete, with small incompletenesses coming primarily from (1) 
galaxies missed because of mechanical spectrograph constraints (~ 7 percent; Blanton et al. 2003), 
which does lead to a slight under-representation of high-density regions, and (2) spectra in which 
the redshift is either incorrect or impossible to determine (< 1 percent). In addition, there are 
some galaxies (~ 1 percent) blotted out by bright Galactic stars, but this incompleteness should 
be uncorrelated with galaxy properties. 

For each galaxy in the sample there is a computed inverse selection volume 1/V max , where 
Vmax is the total comoving volume in which the galaxy could lie and still be included in the sample, 
accounting for the flux limits, surface brightness limits, redshift limits, and completeness as a 
function of angle (as in Blanton et al. 2003a). 

We employ an environment overdensity estimate <5i X 8 based on the SDSS spectroscopic sample; 
it is a measure of the three-dimensional redshift-angle space number density excess around each 



-4 - 



galaxy. The comoving transverse distances and comoving line-of-sight distances (e.g., Hogg 1999) 
are computed between each spectroscopic galaxy and its neighboring spectroscopic galaxies (not 
attempting to correct for peculiar velocities). For each galaxy, neighbors within a cylinder with a 
transverse comoving radius of 1 h^ 1 Mpc and a comoving half-length of 8 h^ 1 Mpc in this space 
are counted; the result is divided by the no-clustering prediction made from the galaxy luminosity 
function (Blanton et al. 2003b) and the sample selection critera, and unity is subtracted to produce 
the overdensity estimate <5i X 8- A galaxy in an environment with the cosmic mean density has 
<5ix8 = 0. The statistical properties of the estimator <5i X 8 vary with redshift because the density of 
the sample varies with redshift, but the sample considered here is selected to be only in the narrow 
redshift interval 0.08 < z < 0.12. It is worthy of note that even galaxies with <5i x g > 50 (the highest 
density subsample used below) are not found exclusively inside the cores of rich clusters. 



3. Results 

Figure 1 shows the distribution in color and magnitude for the entire sample, and subsamples 
chosen to lie in low density regions (S\ x s < 3.0), high density regions (<5i x s > 7.0), and very high 
density regions (<5i X 8 > 50), and subsamples chosen to have low Sersic index (n < 2.0) and high 
Sersic index (n > 2.0). Overplotted on all panels is the same color-magnitude relation, a linear 
fit to the dependence of color on magnitude for the all-environment, high Sersic subsample, with 
sigma-clipping at 2 a in the color direction. The apparent bimodality of the galaxy population 
in Figure 1 and the value of the Sersic index for separating populations has been made evident 
previously (Blanton et al. 2003a). 

The distribution of galaxies is different in the different subsamples. The most striking difference 
is that the bulge-dominated (n > 2.0) galaxies are much redder than the disk-dominated galaxies. 
Another striking difference is that the bulge-dominated galaxies are over-represented in the higher 
density regions. Another is that the very most luminous (Mo.ij ~ —22.8 mag) galaxies are over- 
represented in the higher density regions. 

Among the bulge-dominated (n > 2.0) galaxies it is just as striking that the color-magnitude 
relation is independent of the environmental overdensity. There are environmental differences in 
the distribution of galaxies in the color-magnitude plane, but the most likely color for a bulge- 
dominated galaxy of any given luminosity is not much redder in the highest density environments 
than it is in the lowest. 

Figure 2 shows the color-direction residuals away from the linear fit, for the high Sersic-index 
galaxies in low, high, and very high density regions. The distributions differ in the abundance of 
blue galaxies (and therefore in the shape of the residual distribution), as noted previously (Pimbblet 
et al. 2002; Kuntschner et al. 2002). All three overdensity subsamples show strong red peaks, with 
similar modes and dispersions. In detail, the highest density subsample shows the narrowest and 
reddest peak, with a mode that is slightly redder (by < 0.02 mag in 01 [g — r]) than that of the 
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lowest density subsample. 

This independence of red galaxy colors on environment is consistent with a previous SDSS 
result (Bernardi et al. 2003c) which was based on a restricted sample of early-type galaxies, selected 
on the basis of radial profile, ellipticity, color, and lack of emission lines. Small variations with 
environmental density have been found in the mean strengths of narrow spectral features in red 
galaxies (Eisenstein et al. 2003), but evidently they are not strong enough to affect the broad-band 
colors significantly. 

There are several other trends visible in Figure 1. The bulge-dominated galaxies show a 
larger "tail" of blue outliers in lower density regions; some of the galaxies in this tail are post- 
starburst galaxies (Quintero et al. 2003); some of the others are LINERs (work in preparation). 
The high-density regions contain a small population of red, exponential, sub-L* galaxies (L* is at 
Mo.ij ~ —20.8 mag, Blanton et al. 2003b) not seen at lower densities. This population is visible in 
previous SDSS work (Hogg et al. 2003). 

The slope (color over absolute magnitude) of the color-magnitude relation plotted in Figure 1 
is —0.022 mag mag -1 , and the HWHM of the relation visible in Figure 2 is < 0.05 mag (expressed 
as a limit because calibration errors may contribute significantly). If high luminosity red galaxies 
are made from the mergers of lower luminosity galaxies without further star formation, the color- 
magnitude relation is difficult to explain. On the other hand, the major merger of two equal- 
luminosity galaxies on the color-magnitude relation does not, by itself, create a noticeable outlier 
in this dataset. 

The fact that the distribution of bulge-dominated galaxy luminosities is a strong function 
of environment is emphasized in Figure 3, which shows the luminosity distributions of the high 
Sersic (n > 2.0) galaxies in the three environment subsamples. The abundance of high luminosity 
galaxies increases rapidly with environmental overdensity. This increase in average luminosity with 
increasing overdensity is related to (but not identical with) the previously shown increase in average 
overdensity with increasing luminosity (i.e., the converse relation) for these galaxies (Hogg et al. 
2003). 

A change of ~ 0.01 mag in 0-1 [g — r] corresponds to ~ 0.8 Gyr of aging or ~ 0.04 dex in 
metallicity for an old (~ 10 Gyr), solar metallicity stellar population (e.g., Chariot et al. 2003). 
Of course it is possible that there are compensating effects between metallicity and age (e.g., 
Kuntschner et al. 2002), although no evidence for such a trade-off appears in the analysis of SDSS 
spectra; what appears is evidence that small age effects, small metallicity effects, and small relative 
abundance effects (e.g., alpha enhancement) all come in at comparable levels (Eisenstein et al. 
2003). 

It is difficult to imagine how high and low density regions of the Universe can evolve in their 
different ways without affecting the color-magnitude relation of the bulge-dominated galaxy pop- 
ulation, in color or slope, by more than what is observed here. 
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Fig. 1. — The top left panel shows the distribution in absolute magnitude Mo.i i in the i band 
and 0,1 [g — r] color for the entire sample. The columns show subsamples cut in overdensity 5i x s 
(a mean density environment has <5i X 8 = 0). The rows show subsamples cut in Sersic index n (an 
exponential disk has n = 1 and a de Vaucouleurs profile has n = 4). In each panel, the greyscale 
monotonically represents the abundance of sample galaxies in the two-dimensional space of color 
and magnitude and the contours enclose 52.0, 84.3, and 96.6 percent of the sample. This represents 
properties of the sample; the data have not been weighted by l/V max . Overplotted on all panels 
is the same straight, solid line showing the best-fit color-magnitude relation from the lower-left 
panel, fit as described in the text. 
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Fig. 2. — The distribution of color residuals (observed minus fit) relative to the linear fit plotted in 
Figure 1 for the high sersic samples at low and high density. The histograms have been arbitrarily 
re-normalized to have similar peak heights. 
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Fig. 3. — The 1/Vmax-weighted distribution function of galaxy absolute magnitudes for the high 
Sersic (n > 2.0) galaxies for the different density environments. The error bars show Poisson 
uncertainties only. 



